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SECTION  1 
INTRODUCTION 

The  research  on  this  project  has  been  concerned  with  developing 
techniques  for  the  study  of  liquid  impact,  with  assessing  the  behaviour 
of  various  infrared  transparent  solids  and  most  recently  with  analytic  and 
experimental  work  on  hard  coatings. 

As  discussed  in  detail  in  our  previous  interim  reports  (1,2),  we 
use  three  different  techniques  for  studies  of  high  velocity  liquid/solid 
impact.  The  first  involves  projecting  specimens  of  up  to  25.4  mm  diameter 
at  stationary  drops.  The  second  fires  a  liquid  jet  at  a  stationary  target. 
The  third  uses  two-dimensional  configurations  (discs,  wedges  of  liquid) 
which  are  impacted.  The  first  approach  is  nearest  to  the  practical  situa¬ 
tion  in  liquid  impact,  the  second  has  distinct  advantages  in  its  ease  of 
operation,  low  construction  cost  and  the  velocity  range  which  can  be  covered. 
The  final  approach,  using  the  two-dimensional  configuration,  allows  proces¬ 
ses  occurring  inside  the  impacted  liquid  to  be  followed  by  high-speed  pho¬ 
tography,  and  is  nearer  situations  which  can  be  theoretically  analyzed (3 ,4)  . 

In  order  to  assess  impact  damage  quantitatively,  we  measure  ’residual 
strengths'.  This  involves  impacting  a  specimen  under  known  conditions  and 
then  measuring  the  strength.  A  hydraulic  apparatus,  which  stresses  2" 

(ca.  51  mm)  diameter  disc  specimens,  was  developed  for  this  (5).  The  tech¬ 
nique  has  subsequently  been  modified  and  improved  (6).  The  technique  can 
be  used  for  a  range  of  disc  sizes;  recently  we  have  designed  an  apparatus 
for  stressing  discs  of  1"  (25.4  mm)  diameter.  This  is  proving  useful  where 
the  amount  of  specimen  material  is  limited  or  expensive. 

Our  earlier  work  has  shown  that  it  is  possible  to  relate  drop  and  jet 
impact  in  a  useful  manner,  and  to  use  the  concept  of  an  'equivalent  drop 
size'  for  jet  impacts.  Our  high-speed  camera  work  (1)  showed 

that  our  apparatus  produces  jets  which  are  coherent  and  with  smooth,  curved 
front  faces.  Further,  since  the  high  pressures  of  liquid  impact  occur  in 
the  first  instant  of  impact,  jet  and  drop  impact  can  be  made  'equivalent' 
to  a  reasonable  approximation. 

The  present  report  starts  by  describing  an  analysis  of  'residual' 
strength  curves.  It  is  shown  that  it  is  possible  to  correlate  the  various 
regions  of  the  curves  with  the  response  of  flaws  in  the  material  to  the 
impact  stresses.  The  analysis  predicts  the  shape  of  the  experimental 
curves  very  satisfactorily.  Examples  are  given  of  the  effect  on 
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changes  in  the  flaw  size  distribution  and  changes  in  on  the  material 
response  to  liquid  impact 

Many  of  the  materials  which  have  attractive  IR  transmission  proper¬ 
ties  are  mechanically  weak.  Much  of  our  recent  research  effort  has  been 
concerned  with  assessing  the  possibly  beneficial  effects  of  a  thin 
of  hard  coating.  Sections  3  and  4  describe  the  results  of  finite  element 
analyses  for  the  cases  of  (i)  Hertzian  (ii)  uniform  pressure  loading  of- 
substrates  with  a  hard  coating.  Variables  are  the  elastic  properties  of 
the  coating  and  its  thickness. 
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SECTION  II 

LIQUID  IMPACT  DAMAGE  IN  BRITTLE  SOLIDS 

1.  INTRODUCTION 

In  a  previous  report  (2)  we  have  described  the  damage  in  soda-lime 
glass  due  to  single  and  multiple  jet  impact.  The  damage  was  quantified 
by  measuring  the  post-impact  strength  of  specimens  using  a  specially  des¬ 
igned  hydraulic  bursting  technique  (5,6).  A  typical  'residual  strength' 
curve  for  single  impact  on  a  brittle  material  is  shown  in  Figure  1. 

Three  stages  can  be  distinguished  in  this  curve: 

i)  at  low  impact  velocities  no  damage  can  be  detected  and  the  fracture 
stress  of  the  specimen  is  as  that  of  unimpacted  specimens. 

ii)  at  moderate  impact  velocities  the  average  fracture  stress  decreases 
rapidly  with  increasing  impact  velocity.  In  this  region  there  is  a  bi- 
modal  distribution  in  the  results:  some  specimens  have  a  low  fracture 
stress  while  others,  impacted  at  the  same  velocity,  fail  at  stresses  com¬ 
parable  to  those  for  unimpacted  specimens.  The  percentage  of  specimens 
with  low  fracture  stresses  increases  continually  with  impact  velocity  as 
is  illustrated  in  the  failure  probability  curve  shown  in  Figure  2.  It 
should  be  emphasized  that  earlier  work  by  Rickerby  (7)  using  high-speed 
photography  has  shown  that  the  variability  in  the  results  is  not  caused 
by  jet  instabilities  but  is  due  to  local  variations  in  the  surface  flaw 
size  population.  This  explanation  is  supported  by  a  new  theoretical 
analysis  of  the  impact  damage  which  is  described  later  in  this  section. 

iii)  at  high  impact  velocities  the  average  fracture  stress  still  decreases 
continually  with  impact  velocity  but  at  a  much  lower  rate.  All  specimens 
show  extensive  damage. 

It  has  been  shown  by  Bowden  and  Field  (8)  that  the  jet/drop  impact 
damage  in  brittle  materials  is  due  to  the  interaction  of  the  Rayleigh 
surface  waves  with  pre-existing  surface  defects  and  the  ensuing  extension 
thereof.  In  the  previous  report  (2)  we  have  calculated  the  growth  of  sur¬ 
face  flaws  due  to  stress  pulses  simulating  the  Rayleigh  wave  generated  by 
the  impact  as  a  first  attempt  to  calculate  the  velocity  dependence  of  jet/ 
drop  impact  damage. 

In  this  report  we  have  extended  and  improved  the  previous  model  sig¬ 
nificantly  by  taking  into  account  the  statistical  nature  of  the  flaws  in 
brittle  materials. 
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2.  A  MODEL  FOR  LIQUID  IMPACT  DAMAGE 

We  now  consider  in  some  detail  the  various  parameters  which  deter¬ 
mine  the  impact  stresses  and  the  resulting  damage  and  combine  them  to  give 
a  simple  quantitative  model  for  liquid  impact  damage. 

a.  Stress  Wave  Intensity 

Swain  and  Hagan  (9)  have  attempted  to  measure  the  magnitude  and 
time  dependence  of  the  Rayleigh  surface  wave  produced  by  an  impacting  jet 
using  piezo-electric  crystals.  They  were  unable  to  determine  the  absolute 
magnitude  of  the  pulse,  but  showed  that  the  Rayleigh  wave  can  be  approxi¬ 
mated,  to  a  reasonable  accuracy,  by  a  triangular  pulse.  As  a  first  ap¬ 
proximation  they  used  the  following  expression  for  the  amplitude  of  the 
Rayleigh  wave 

a  =  a  Z  V.  (1) 

max  1 

where  a  is  a  constant,  Z  the  acoustic  impedance  of  the  liquid  and  the 
impact  velocity.  The  above  equation  was  derived  by  drawing  an  analogy 
with  the  case  of  quasi-static  Hertzian  contact  where  the  maxium  radial 
tensile  stress  depends  linearly  on  the  average  pressure  over  the  contact 
zone. 


b.  Dynamic  Crack  Growth 

In  the  case  of  a  dynamically  loaded  crack  the  stress  intensity  factor 
is  initially  a  function  of  time  (10-12).  For  a  stationary  crack  the 
diffraction  of  a  stress  wave  by  a  crack  of  finite  dimensions  leads  to  a 
damped  oscillatory  behaviour  of  the  stress  intensity  factor  around  the 
quasi-static  value. 

The  stress  intensity  factor  for  a  non-stationary  crack  has  been 
analysed  in  several  studies  (13-18).  The  equation  for  the  stress  intensity 
factor  derived  by  Eshelby  (14)  is  particularly  applicable  to  the  problem 
of  propagating  small  flaws  and  is  given  by 

=  m  a(7ra)^(l  -  V/Cr)^  U  (2) 


where  m  is  a  geometrical  constant,  a  the  initial  flaw  size,  V  the  current 
crack  velocity,  Cr  the  Rayleigh  wave  velocity  and  U  is  given  by  the 
following  series 


U  = 


1  + 


8 


g  -  a  _  iA  fC  -  a, 
a  26  '  a  ^ 


(3) 


where  c  is  the  current  crack  size.  Evaluation  of  this  series  shows  that 
for  c  ^  20a  U  can  be  accurately  approximated  (16)  to 
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u 


(4) 


The  current  crack  velocity  in  the  above  equation  can  be  derived  from  the 
experimentally  determined  dependence  of  the  crack  velocity  on  the  dynamical 
stress  intensity  factory  (19) 


V 


V  (1 
max 


(5) 


where  V  is  the  maximum  crack  velocity, 
max  ^ 

c.  Flaw  Statistics 

The  variability  of  the  strength  of  nominally  identical  brittle 
specimens  is  due  to  a  distribution  in  the  size  of  the  'Griffith'  flaws. 
Hunt  and  McCartney  (20)  have  derived  an  expression  relating  the  failure 
probability  to  the  flaw  size  distribution.  A  simplified  probability  den¬ 
sity  function  for  the  flaw  dimensions  has  been  derived  by  Jaytilaka  and 
Trustrum  (21)  from  experimental  data  and  is  given  by 


f  (c) 


-  n-1  , 

fS _  c'” 

(n-2)! 


(6) 


where  c^  is  the  size  of  the  most  probable  flaw  dimension.  The  parameter 
n  is  related  to  the  Weibull  parameter  m  (22). 


d.  The  Model 


In  order  to  calculate  the  residual  strength  curve  for  soda-lime 
glass  the  following  procedure  has  been  followed. 

i)  In  a  circular  specimen  flaw  positions  were  generated  at  random.  The 
total  number  of  flaws  per  specimen  was  between  20  and  40. 

ii)  A  flaw  dimension  was  allocated  to  each  of  the  flaw  positions.  All 
flaws  were  taken  as  normal  to  the  radial  vector  from  the  centre  of  the 


specimen.  The  flaw  size  distribution  used  in  the  calculation  is  plotted 
in  Figure  3.  The  most  frequent  crack  dimension  was  taken  as  12  pm;  99% 
of  the  flaws  have  dimensions  smaller  than  50  pm. 


iii)  The  maximum  flaw  size  in  the  specimen  was  used  to  calculate  the  pre¬ 
impact  strength  of  this  specimen. 

iv)  The  impact  was  simulated  by  a  stress  pulse,  i.e.  the  Rayleigh  wave, 

moving  radially  outwards  from  the  centre  of  the  specimen.  In  order  to 
simplify  the  calculations  the  stress  pulse  was  taken  as  rectangular.  The 
amplitude  of  this  pulse  decayed  as  R  i.e.  the  same  decay  rate  as 

the  Rayleigh  wave,  when  moving  outwards  over  the  specimen.  The  initial 
duration  of  the  stress  pulse  was  taken  as  0.1  ps.  The  duration  of  the  pulse 


5 


1 


was  increased  linearly  with  radial  distance. 

(v)  At  each  flaw  position  the  dynamical  stress  intensity  factor  was  cal¬ 
culated  according  to  equation  2.  For  flaws  where  exceeeded  crack 
growth  was  calculated  using  an  interation  procedure  in  which  the  stress 
pulse  was  divided  in  100  intervals  of  equal  duration.  For  each  interval 
the  current  crack  velocity  was  calculated  using  equation  5.  From  this 
crack  velocity  the  crack  growth  during  the  interval  was  calculated  and 
added  to  the  current  crack  dimension  calculated  in  the  previous  interval. 

In  the  calculations  the  appropriate  values  for  K  (0.75  MNm“3/2)  and 
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maximum  crack  velocity  (1500  m  s~*)  have  been  used. 

(vi)  From  the  maximum  flaw  size  after  impact  the  post-impact  strength  of 
this  specimen  was  calculated.  By  comparing  the  pre  and  post-impact  strength 
of  the  specimen>the  occurence  of  impact  damage  could  be  determined. 

(vii)  The  above  procedure  was  repeated  for  a  100  specimens  under  identical 
impact  conditions.  For  each  impact  velocity  the  average  fracture  stress 
and  the  probability  that  the  fracture  stress  is  reduced  by  the  impact 
were  calculated. 

(viii)Finally  the  initial  stress  pulse  amplitude,  i.e.  the  impact  velocity 
(equation  1),  was  stepwise  increased  and  the  whole  procedure  was  repeated 
for  each  pulse  amplitude. 

The  results  of  the  calculations  have  been  plotted  in  the  form  of  a 
residual  strength  curve  in  Figure  4.  A  very  good  qualitative  agreement 
with  the  residual  strength  curve  for  single  impact  is  obtained.  In  parti¬ 
cular  the  variation  in  the  amount  of  scatter  in  the  fracture  stress  as  a 
function  of  the  impact  velocity  is  well -reproduced  in  the  analysis.  From 
the  analysis  a  failure  probability  curve  can  be  obtained  which  is  also 
shown  in  Figure  4  (broken  line).  Once  again  a  good  qualitative  agreement 
with  the  measured  probability  curve  is  obtained.  The  model  has  also  been 
used  to  calculate  the  effect  of  multiple  impacts  on  the  resulting  damage 
and  the  same  conclusions  as  from  the  experimental  results  can  be  made, 
namely  an  almost  constant  threshold  velocity,  a  reduction  of  the  width 
of  the  transition  region  and  a  continuous  increase  in  the  post-impact  flaw 
size  at  high  impact  velocities. 

The  model  can  also  be  used  to  illustrate  the  effect  of  abrasion  on 
the  residual  strength  curve.  For  this  purpose  a  new  set  of  specimens  with 
twice  the  density  of  flaws  was  generated.  Furthermore,  the  size  of  the 
most  frequent  defect  was  increased  to  25  pm  with  99%  of  the  defects  smaller 
than  too  pm.  The  impact  parameters  and  material  properties  were  not 
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changed.  The  residual  strength  curves  for  the  'as  received'  (open  circles) 
and  the  'abraded'  specimens  (closed  circles)  are  plotted  in  Figure  5.  Due 
to  the  abrasion  not  only  the  initial  strength  of  the  specimens  is  reduced 
but  also  the  threshold  velocity  and  the  post  impact  strength  at  high  impact 
velocities.  The  failure  probability  curve  (not  shown)  shows  a  much  smal¬ 
ler  transition  region  for  the  'abraded*  specimens.  The  above  predictions 
are  in  very  good  agreement  with  experimental  data  for  jet  impact  studies 
(23).  Furthermore,  Wiederhorn  and  Lawn  (24)  noted  similar  effects  for 
solid  particle  impacts  on  abraded  and  unabraded  glass  specimens. 

Finally,  in  Figure  6  the  beneficial  effect  of  increased  toughness  of 

the  specimens  is  illustrated.  The  solid  line  indicates  the  residual  strength 
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curve  for  a  material  with  K^_  =  1.5  MKn  (i.e. twice  that  of  the  reference 
material).  The  other  impact  and  material  parameters  are  as  for  the  refer¬ 
ence  conditions.  It  should  be  mentioned  that  such  an  increase  in  impact 
resistance  can  also  be  obtained  by  polishing  the  specimens,  i.e.  producing 
smaller  flaws  with  a  lower  density.  However,  in  practice  polishing  is  not 
a  feasible  technique  for  improving  the  rain  erosion  resistance  because  of 
the  problems  associated  with  maintaining  the  surface  finish. 


SECTION  III 

THE  EFFECT  OF  THIN  HARD  COATINGS  ON  THE  HERTZIAN  STRESS  FIELD 


I-  INTRODUCTION 

The  elastic  contact  between  brittle  solids  and  spherical  indenting 
or  impacting  particles  may  lead  to  the  formation  of  ring  and  Hertzian 
cone  cracks.  These  cracks,  which  are  due  to  radial  tensile  stresses  in  a 
shallow  surface  layer  around  the  contact  area,  cause  surface  degradation 
and  can  lead  to  significant  losses  in  strength  and  transparency.  This 
form  of  contact  damage  has  restricted  the  use  of  brittle  materials  in  many 
high  performance/reliability  applications.  However,  since  ceramics  and 
glasses  possess  a  wide  range  of  useful  properties  which  are  not  (or  not 
easily)  obtainable  from  more  ductile  materials,  the  integration  of  brittle 
components  in  engineering  structures  is  now  frequently  desirable  if  not 
essentia] .  One  way  of  reducing  the  damage  due  to  static  or  dynamic  con¬ 
tact  is  by  protecting  the  surface  with  a  coating.  The  application  of  com¬ 
pliant  protective  coatings  is  well-established  (e.g.  plastic-coated  bottles 
and  neoprene-coated  aircraft  radomes)and  has  been  examined  in  detail  Pitts (25) ; 
Dhaliwahl  and  Rau(26);  Conway,  Lee  and  Bayer  (27);  Conway,  Engel  and  Lee 
(28);  Chen  and  Engel  (29);  Matthewson  (30,31).  Of  particular  interest  is 
recent  theoretical  work  by  Matthewson  (30,31)  in  which  the  indentation 
stress  field  is  calculated  for  a  compliant  coating  bonded  to  a  rigid  sub¬ 
strate.  It  is  shown  that  the  presence  of  a  coating  of  appropriate  thick¬ 
ness  can  significantly  reduce  the  maximum  radial  tensile  stress  in  the  sub¬ 
strate  and  that  incompressible  materials  are  most  effective.  Much  of  the 
protection  is  however  lost  when  coating  debonding  occurs. 

Although  metallic  and  elastomeric  coatings  have  proved  their  useful¬ 
ness  there  are  now  several  practical  situations  (e.g.  infrared  transparent 
aircraft  windows)  where  these  traditional  coating  materials  cannot  be  used 
because  of  unsuitable  electromagnetic  properties.  In  order  to  determine 
the  feasibility  of  protecting  brittle  substrates  by  a  thin  hard  coating  and 
to  direct  the  development  of  new  ceramic  coating  materials,  a  static  finite 
element  model  (using  the  Package  for  Automatic  Finite  Element  Calculations 
described  by  Henshell  (32))  is  developed  in  which  an  undeformable  sphere 
is  loaded  onto  a  coated  substrate.  In  the  model  both  the  substrate  and 
coating  stresses  are  calculated  as  a  function  of  the  elastic  properties  of 
the  coating  material  and  the  coating  thickness  and  are  compared  with  the 
analytical  solutions  for  the  indentation  stress  field  in  a  uniform  halfspace. 
The  model  is  described  in  Section  3.3. 


2. 


THE  HERTZIAN  STRESS  FIELD 

The  analytical  solution  of  the  elastic  stress  field  generated  by  two 
spherical  bodies  in  frictionless  contact  has  been  determined  in  general 
terms  by  Hertz  (33)  and  has  further  been  developed  by  Huber  (34),  Fuchs 
(35),  Huber  and  Fuchs (36)  and  Morton  and  Close  (37).  Later  analyses  have 
considered  the  effect  of  friction,  Spence  (38),  Johnson,  O'Connor  and 
Woodward  (39), and  surface  roughness.  Greenwood  and  Trip  (40),  on  the  stress 
field  generated.  In  the  present  work  the  results  for  the  indentation  stress 
field  in  composite  structures  are  compared  with  the  analytical  solutions 
for  the  radial,  axial  and  circumferential  stresses  as  derived  by  Huber  (34). 

For  a  Hertzian  contact  condition,  in  which  a  ball  of  radius  R  is 
pressed  onto  a  flat  uniform  half space  by  a  force  F,  the  radius,  a,  of  the 
contact  zone  is  given  by  the  well“known  equation 
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where  E|,  and  are  the  Young's  modulus  and  the  Poisson's  ratio 

of  the  halfspace  and  the  indenter  respectively.  The  semi-ellipsoidal  pres¬ 
sure  distribution  over  the  surface  of  the  contact  zone  is  given  by 
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The  tangential  displacements  of  the  surface  due  to  this  pressure  distribu¬ 
tion  are  radially  inwards  and  lead  to  a  region  of  tensile  radial  stresses 
around  the  edge  of  the  contact  zone.  The  radial  stress  distribution,  both 
within  and  outside  the  contact  zone,  is  given  by 


3F  1 

1  -  2v  a^ 

i-^-V 

♦  f- 

2  * 
2TTa^  j 

■>  2 

3  r 

1 

2  2 
)  z 


z 

7u' 


(1  -  V  )  u  /u  a 

+  (1  +  v)  —  arctg  —  -  2 

a  /u 


2  ^ 

a  +  u 


(9) 


where  r  is  the  radial  coordinate,  and  u  the  positive  root  of  the  quad¬ 
ratic  equation  :  u^  +  a^u  -  r^  -  z^  =  0,  At  the  surface  and  within 
the  contact  zone  the  radial  stress  distribution  is  given  by 
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The  radial  stress  is  compressive  within  the  contact  zone  but  is  tensile 
within  a  shallow  surface  layer  outside  the  contact  zone.  In  the  case  of 
elastic  ball  indentations  on  brittle  materials,  the  tensile  radial  stress 
outside  the  contact  zone  is  the  most  important  of  the  stresses  since  it 
can  lead  to  the  propagation  of  surface  flaws  which  cause  strength  reduc¬ 
tions,  material  removal  and  surface  degradation. 

The  special  dependence  of  the  circumferential  and  axial  stresses 
are  given  by 


(ID 


(12) 


where  the  parameters  have  the  same  meaning  as  in  the  previous  equations. 
Note  that  equation  1 1  is  the  corrected  form  of  the  original  Huber  equation 
(41). 

3.3  FINITE  ELEMENT  MODEL  FOR  SPHERICAL  INDENTATIONS  ON  COATED  SURFACES 

In  this  analysis  for  the  static  indentation  of  a  coated  substrate  by 
an  undeformable  spherical  indenter  the  effect  of  the  thickness  and  the 
elastic  properties  of  the  coating  are  investigated.  Throughout  the  analy¬ 
sis  the  properties  of  the  substrate  are  held  constant  and  equal  to  the  pro¬ 
perties  of  glass  as  given  in  the  PAFEC  materials  module  namely  :  Young's 
modulus  61  GPa  and  Poisson's  ratio  0.25.  The  effect  of  the  coating  thick¬ 
ness  has  been  investigated  for  thicknesses  up  to  20%  of  the  contact  radius. 
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i.e.  only  thin  coatings  have  been  modelled.  The  Young's  modulus  of  the 
coating  material  has  been  varied  between  I  and  10  times  that  of  the  sub¬ 
strate  material.  This  range  covers  nearly  all  possible  hard  coating  mater¬ 
ials  (e.g.  the  Young's  modulus  of  diamond  is  about  15  times  as  high  as  that 
of  the  chosen  substrate).  The  Poisson's  ratio  of  the  coating  has  been  var¬ 
ied  between  0.20  and  0.30.  Once  again,  this  range  covers  the  majority  of 
the  potential  hard  coating  materials. 

The  finite  element  distribution  used  is  shown  in  Figure  7a.  The 
elements  are  divided  into  two  groups  :  one  with  the  (fixed)  substrate 
properties  and  one  with  the  (variable)  coating  properties.  The  elements 
with  the  coating  properties  are  shaded  in  Figure  7b.  The  thickness,  d, 
of  the  coating  is  very  small  compared  with  the  dimensions  of  the  total 
structure  (less  than  2%),  so  that  boundary  effects  do  not  interfere  with 
the  stress  field  calculated.  It  has  been  verified  that  in  the  case  of  a 
uniform  half space  (i.e.  the  coating  elements  have  the  same  elastic  proper¬ 
ties  as  the  substrate)  the  finite  element  results  for  the  radial,  axial  and 
circumferential  stress  agree  with  the  appropriate  analytical  solutions. 

For  the  chosen  element  distribution,  the  difference  between  the  two  solu¬ 
tions  is  always  less  than  0.5%. 

In  the  analysis  only  vertical  displacements,  w,  are  prescribed  to 
the  nodes  within  the  chosen  contact  zone  according  to 

2  2 

w  =  —  -  —  (13) 
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In  order  to  model  the  case  of  elastic  contact  accurately,  the  ratio  of 
the  contact  radius  to  the  indenter  radius  is  chosen  as  1  to  100.  Because 
of  the  discontinuity  in  the  radial  and  circumferential  stress  across  the 
interface  the  reported  stresses  at  the  interface  are  those  calculated  for 
the  interface  elements  with  substrate  properties.  Such  an  approach  gives 
a  true  indication  of  the  effect  of  the  various  coating  parameters  on  the 
substrate  stress  field. 

In  the  following  sections  the  calculated  stresses  are  normalised  to 
appropriate  reference  stresses.  The  normalisation  was  made  such  that 
positive  values  indicate  tensile  stresses  and  negative  values  compressive 
stresses . 

3.4  RESULTS  FOR  AN  EXAMPLE  PROGRAM 

Before  considering  the  main  results,  the  results  of  a  particular 


PAFEC  program  are  shown  and  discussed.  In  this  example  program  =  610  GPa, 
=  0.20,  E^  =  61  GPa,  =  0.25  and  the  ratio  of  the  coating  thickness  to 
the  contact  radius  =  0.05.  The  subscripts  c  and  s  refer  to  the  coating  and 
substrate  material  respectively. 

The  calculated  axial  stress  distribution  at  the  interface  is  shown 
in  Figure  8.  The  axial  stress  is  normalised  to  the  maximum  axial  stress  Pq 
calculated  analytically  for  an  indentation  of  the  same  contact  radius  on  a 
uniform  half space  with  the  substrates  properties.  The  finite  element  res¬ 
ults  are  indicated  in  this  figure  by  the  closed  circles.  The  dashed  line 
gives  the  analytical  solution  for  the  radial  stresses  in  a  halfspace  with 
substrate  properties  at  the  depth  of  the  interface.  The  figure  shows  that 
for  the  same  contact  radius  the  axial  stresses  are  larger  due  to  the 
presence  of  the  stiff  coating,  as  would  be  expected  intuitively.  The 
maximum  axial  stress  is  about  10%  higher  than  for  Hertzian  contact  on  a 
uniform  half space.  The  spatial  distribution  of  the  axial  stress  in  the 
contact  zone  is  however  not  changed  significantly  by  the  coating.  This  is 
illustrated  by  the  good  fit  between  the  finite  element  results  and  the  solid 
line,  which  is  given  by 
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where  P  is  the  maximum  axial  stress  on  the  coating  surface  as  calculated 

by  the  finite  element  program.  The  equation  above  is  identical  to  the 

equation  for  the  axial  stress  (equation  8)  derived  by  Huber  (34),  apart 

* 

from  the  scaling  factor  P  . 

The  circumferential  stress  at  the  interface  is  shown  in  Figure  9.  The 
same  normalisation  procedure  as  for  the  axial  stress  field  is  used.  The 
finite  element  results  are  indicated  by  closed  circles.  The  dashed  line 
shows  the  Hertzian  solution  for  the  hoop  stress  at  the  depth  of  the  inter¬ 
face.  The  solid  line  is  identical  to  the  Hertzian  equation  for  the  circum¬ 
ferential  stress  (equation  11),  apart  from  the  scaling  factor  P  .  This 
figure  shows  that  the  coating  reduces  the  circumferential  stresses  consider¬ 
ably  both  within  and  outside  the  contact  zone.  The  reduction  is  due  to  the 
restricted  displacements  of  the  substrate  surface  because  of  the  higher 
stiffness  of  the  rigidly  adhering  coating. 

The  radial  stress  at  the  interface  is  shown  in  Figure  10.  The  finite 
element  results  are  again  indicated  by  closed  circles.  The  dashed  line  is 
the  analytical  solution  for  the  radial  stress  in  a  uniform  halfspace  at  the 
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depth  of  the  interface.  The  solid  line  is  a  similar  modification  of  the 
original  Huber  solution  as  has  been  made  for  the  axial  and  circumferential 
stresses.  The  figure  shows  clearly  that  the  maximum  radial  stress  in 
the  substrate  is  considerably  reduced.  This  indicates  that  hard  coatings 
can  be  used  to  protect  the  substrate  against  damage  due  to  localised  load¬ 
ing.  The  reduction  of  the  radial  stresses  can  be  explained  in  the  same 
way  as  the  reduction  in  the  circumferential  stresses,  namely,  by  a  reduc¬ 
tion  of  the  radial  displacements  at  the  interface.  This  reduction  in 
radial  displacements  is  shown  in  Figure  1 1 ,  which  shows  the  radial  dis¬ 
placements  (normalised  to  the  maximum  penetration  of  the  indenter)  at  the 
interface  for  the  composite  structure  and  the  radial  displacements  at  the 
depth  of  the  interface  for  a  uniform  halfspace  (solid  line). 

Figure  12  shows  the  interfacial  shear  stress  in  a  radial  direction 
both  for  the  composite  structure  and  for  a  'substrate  halfspace'  (solid 
line).  For  the  present  coating  parameters,  the  maximum  shear  stress  is 
more  than  doubled.  This  is  important  since  interfacial  shear  stresses  of 
significant  magnitude  can  lead  to  the  coating  debonding,  in  which  case  the 
protection  provided  by  the  coating  is  greatly  reduced. 

Examination  of  the  stresses  at  the  coating  surface  has  shown  that  the 
axial  stress  or  normal  pressure  at  the  coating  surface  is  nearly  identical 
to  that  at  the  interface  as  can  be  expected  for  any  thin  coating.  The  max¬ 
imum  pressure  at  the  centre  of  the  contact  is  10%  higher  than  for  an  uncoated 
substrate. 

The  effect  of  the  properties  of  the  coating  material  on  the  stresses 
in  the  plane  of  the  coating,  e.g.  the  radial  and  the  circumferential  stress, 
is  however,  very  large.  This  is  illustrated  in  Figure  13  in  which  the  ra¬ 
dial  stress  at  the  coating  surface  has  been  plotted.  The  figure  also  in¬ 
cludes  the  radial  stress  distribution  for  indentations  with  the  same  contact 
radius  on  a  halfspace  with  substrate  and  coating  properties  respectively. 

It  is  shown  that  within  the  contact  zone  the  compressive  stresses  are  much 
higher  than  for  a  substrate  halfspace  but  not  as  high  as  for  a  coating 
half space.  Outside  the  contact  zone  however  the  radial  tensile  stresses 
are  more  comparable  to  those  for  a  coating  halfspace.  The  maximum  tensile 
stress  for  the  present  coating  parameters  is  nearly  equal  to  the  maximum 
axial  stress.  For  the  substrate  half space  the  maximum  tensile  stress  is 
only  16.7%  of  the  central  pressure,  while  for  the  coating  material  this 
ratio  is  20%. 
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3.5  EFFECT  OF  THE  COATING  PARAMETERS  ON  THE  AXIAL  STRESS 


In  Figure  14  the  maximum  pressure,  P  ,  normalised  to  the  maximum 
pressure,  Pq,  for  an  indentation  of  the  same  dimensions  on  the  uncoated 
substrate,  is  plotted  as  a  function  of  the  coating  parameters,  i.e.  the 
coating  thickness,  its  Young's  modulus  and  its  Poisson's  ratio.  Since 
the  pressure  distribution  on  the  coating  surface  is  semi-ellipsoidal, 
this  figure  is  also  a  plot  of  the  normalised  load  versus  the  coating  para¬ 
meters.  The  figure  shows  that  the  maximum  axial  stress  depends  primarily 
on  the  coating  thickness  and  the  Young's  modulus  of  the  coating  material, 
but  generally  the  effect  is  relatively  small  for  these  thin  coatings.  In 
the  figure  the  maximum  pressure  for  an  infinitely  thick  coating  with 
^  =  0.25  (dashed  line)  is  also  included.  In  this  case  the  stress  field 
is  identical  to  that  for  an  indentation  on  a  uniform  half space.  There¬ 
fore,  for  large  coating  thicknesses  the  maximum  pressure  increases  linearly 
with  the  Young's  modulus  of  the  material. 

The  data  plotted  in  this  figure  can  be  easily  converted  to  show  the 
dependence  of  the  contact  radius  on  the  applied  load  for  these  coated 
substrates.  The  results  are  shown  in  Figure  15.  This  figure  also  shows  the 
analytical  solutions  for  the  load  dependence  of  the  contact  radius  for  an 
uncoated  substrate  and  for  two  coating  halfspaces.  For  the  composite  struc¬ 
tures  investigated,  the  analytical  relation  between  load  and  contact  radius 
given  by  equation  7  holds  reasonably  true,  provided  the  coating  thickness 
is  less  than  10%  of  the  contact  radius. 

3.6  EFFECT  OF  COATING  PARAMETERS  ON  THE  RADIAL  STRESS  IN  THE  SUBSTRATE 

In  the  analysis  presented  so  far  the  stress  fields  in  coated  and  un¬ 
coated  substrates  have  been  compared  for  indentations  with  the  same  contact 
radius.  However,  in  such  a  comparison  the  total  applied  load  on  the  inden- 
ter  varies  with  the  coating  parameters.  In  order  to  study  the  effect  of 
the  coating  parameters  for  a  fixed  load,  all  radial  stresses  are  normalised 
to  the  maximum  pressure  for  that  particular  structure. 

The  radial  stress  distribution  at  the  interface  is  shown  in  Figure  16 
for  four  different  Young's  moduli  of  the  coating  material.  The  ratios  of 
the  Young's  modulus  of  the  coating  to  that  of  the  substrate  are  2,3,5  and  10 
respectively.  The  Poisson's  ratios  of  the  coating  and  substrate  material  and 
the  coating  thickness  are  held  constant  at  0.25  and  O.Ia  respectively.  The 
figure  shows  that  within  the  contact  zone  the  magnitude  of  the  radial  stress 
decre  ises  with  increasing  Young's  modulus  of  the  coating  material.  The 
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coating  also  affects  the  maximum  radial  stress  at  the  interface.  The  re¬ 
ductions  of  the  tensile  stresses  outside  the  contact  radius  are  shown  in 
more  detail  in  Figure  17  for  the  same  set  of  coating  conditions.  For  the 
stiffest  coating  investigated,  the  tensile  component  of  the  radial  stress 
field  has  almost  vanished. 

The  coating  not  only  reduces  the  magnitude  of  the  maximum  radial 
stress  but  also  affects  its  position:  it  is  found  at  increasing  radial 
distance.  Since  crack  propagation  in  brittle  materials  depends  on  both 
the  magnitude  of  the  stress  field  and  the  presence  of  sufficiently  large 
flaws,  the  above  increase  in  radial  distance  increases  the  probability  of 
stressing  a  major  flaw.  However,  since  the  radial  tensile  stress  is 
everywhere  smaller  than  that  for  an  uncoated  substrate,  these  hard  coat¬ 
ings  do  indeed  protect  the  surface. 

The  effect  of  the  various  coating  parameters  on  the  maximum  radial 
stress  obtained  in  the  substrate  is  summarized  in  Figure  18.  It  should  be 
emphasized  that  the  maximum  tensile  stress  in  an  uncoated  substrate  decrea¬ 
ses  very  rapidly  with  depth  below  the  surface.  For  a  correct  interpretation 
of  the  effect  of  the  coating,  the  calculated  maximum  tensile  stress  should 
therefore  be  compared  with  the  maximum  tensile  stress  in  a  uniform  half¬ 
space  at  a  depth  identical  to  the  coating  thickness  (dashed  line) .  The 
reduction  is  most  sensitive  to  the  Young's  modulus  and  thickness  of  the 
coating.  For  relatively  thick  coatings  the  Poisson's  ratio  of  the  coating 
material  is  also  important.  The  figure  shows  that  the  effect  of  the  coating 
thickness  on  the  reduction  of  the  radial  stress  is  largest  for  very  thin 
coatings . 

3.7  EFFECT  OF  COATING  PARAMETERS  ON  THE  RADIAL  STRESS  ON  THE  SURFACE 

In  contact  problems  on  brittle  coating  and  substrate  materials  the 
stresses  at  the  surface  of  the  coating  are  as  important  as  the  stresses 
in  the  substrate.  Failure  of  the  total  structure  may  be  caused  by  the 
propagation  of  flaws  located  either  at  the  surface  or  at  the  interface. 

The  effect  of  the  parameters  on  the  maximum  tensile  stress  at  the  coating 
surface,  which  is  also  the  maximum  tensile  stress  in  the  total  structure, 
is  summarised  in  Figure  19.  These  maximum  tensile  stresses  have  been 
determined  by  extrapolation  of  the  radial  stresses  outside  the  contact 
radius  towards  the  edge  of  the  contact  since  the  program  deals  with  the 
stress  singularity  in  an  imperfect  manner.  The  maximum  stress  increases 
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initially  rapidly  with  coating  thickness  for  all  coatings  analysed.  Unlike 
the  case  of  uniform  half spaces,  the  maximum  tensile  stress  can  now  be  larger 
than  the  maximum  pressure  applied.  At  larger  coating  thicknesses  the  max¬ 
imum  radial  stress  reaches  a  maximum  value  and  then  decreases  again.  This 
behaviour  can  be  explained  as  follows  :  for  thin  coatings  the  effect  of  the 
coating  on  the  stress  normal  to  the  plane  of  the  coating,  i.e.  the  axial 
stress,  is  rather  small.  However,  the  coating  has  already  a  large  influence 
on  the  stresses  in  the  plane  of  the  coating,  i.e.  on  the  radial  and  circum¬ 
ferential  stresses.  The  maximum  possible  stress  at  the  coating  surface  is 
that  for  an  indentation  on  a  coating  half space  with  the  same  contact  size. 
The  assumption  that  the  effect  of  the  coating  thickness  on  the  maximum  pre¬ 
ssure  can  be  ignored,  leads  to  the  following  expression  for  the  maximum 
radial  stress  on  the  coating  surface  (for  a  coating  material  with  a  Pois¬ 
son's  ratio  of  0.25) 

0  =  0. 167  P  E  /E  (15) 

max  o  c  s 

where  P^  is  the  maximum  pressure  for  a  similar  indentation  on  the  uncoated 
substrate.  When  the  coating  has  a  different  Poisson's  ratio  the  argument 
is  slightly  more  complicated,  since  the  effect  of  the  Poisson's  ratio  on 
the  contact  radius  has  to  be  taken  into  account.  From  the  Hertzian  solu¬ 
tion  the  following  threshold  values  are  derived 

a  =  0. 195  P  E  /E  (16) 

max  o  c  s 

for  a  coating  material  with  Poisson's  ration  0.20  and 

0  =  0.  137  P  E  /i;  (17) 

max  o  c  s 

for  a  coating  material  with  Poisson's  ratio  0.30.  These  threshold 
values  are  indicated  in  the  figure  (dashed  lines).  A  good  agreement  be¬ 
tween  the  calculated  maximum  value  of  the  radial  stress  and  the  predicted 
thresholds  is  obtained  for  coatings  with  relatively  low  Young's  moduli. 

In  the  case  of  stiffer  coating  materials  the  assumption  made  is  no  longer 
valid  and  the  predicted  threshold  value  should  be  well  above  the  calculated 
maximum,  as  is  indeed  observed.  For  very  thick  coatings  the  response  of 

the  composite  structure  to  the  imposed  indentation  approaches  that  of  a 

* 

coating  halfspace.  In  that  case  the  ratio  between  given 

by  the  appropriate  Hertzian  solution. 
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3.8  COATING  DEBONDING 

In  the  work  presented  so  far  attention  has  been  focussed  on  the  : 

effect  of  the  coating  properties  on  the  tensile  radial  stress  component 
"  i.e.  on  the  brittle  failure  mode  of  the  composite  structure.  It  is  shown 

'  that  a  stiff  coating  will  reduce  the  radial  displacements  at  the  interface 

and  hence  the  radial  stresses  in  the  substrate.  Since  the  protective 

influence  of  the  coating  is  based  upon  good  adhesion  across  the  interface,  i - 

hardly  any  protection  will  be  offered  by  these  thin  coatings  when  debond-  |. 

ing  occurs.  The  debonding  is  determined  by  the  shear  strength  of  the  | 

interface  and  the  shear  stresses  due  to  the  applied  load.  In  the  case  of  | 

polymer  coatings  it  has  been  observed  that  the  shear  strength  of  the  coat-  ‘ 

ing  also  depends  on  the  magnitude  of  the  normal  component  of  the  load  (42).  ii 

No  such  information  is  available  for  ceramic  coatings  but  the  effect  will  ' 

probably  be  negligible  for  this  type  of  material.  Finally,  the  critical 
conditions  for  coating  debonding  depend  also  on  the  magnitude  of  the  pre¬ 
existing  stresses  at  the  interface.  This  is  an  important  problem  for 
hard  coatings,  in  particular  for  coatings  produced  by  chemical  vapour 
deposition  processes  where  the  deposition  takes  place  at  elevated  tempera¬ 
tures.  On  cooling  down, residual  stresses  will  develop  due  to  differences 
in  the  thermal  expansion  coefficients  of  the  coating  and  substrate  mater¬ 
ials.  The  magnitude  of  these  stresses  increases  with  coating  thickness 
and  therefore  limits  the  maximum  coating  thickness  which  can  be  deposited 
before  spontaneous  coating  debonding  occurs. 

The  interfacial  shear  stress  is  plotted  in  Figure  20  as  a  function 
of  the  coating  thickness  for  constant  coating  properties  (E^  =  10  E^;  1 

v^  =  0.20).  The  maximum  shear  stress  is  only  weakly  dependent  upon  the 
!  coating  thickness  but  seems  to  reach  a  maximum  for  a  coating  thickness 

I  of  about  10%  of  the  contact  radius.  For  very  thin  coatings  the  region 

t 

I  of  high  shear  stress  is  localised  around  the  contact  edge.  For  thicker 

i  coatings  high  shear  stresses  are  generated  over  a  much  larger  agrea. 

I  Qualitatively  the  same  observations  are  made  for  coatings  with  different 

elastic  properties. 

Figure  21  shows  the  calculated  interfacial  shear  stresses  for  three 
coating  stiffnesses  (2,5  and  10  Eg).  The  other  coating  parameters, 
thickness  (0.1a)  and  Poisson's  ratio  (v  =  0.25),  are  kept  constant.  The 
I  maximum  shear  stress  increases  with  the  Young's  modulus  of  the  coating 

I  material.  Furthermore,  high  shear  stresses  are  found  over  larger  areas 

I  for  higher  coating  stiffnesses. 
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Examination  of  the  effect  of  the  Poisson's  ratio  has  shown  that  the 


shear  stresses  are  virtually  independent  of  the  Poisson's  ratio  for  the 
limited  range  investigated;  this  applies  in  particular  to  the  very  rigid 
coatings.  The  same  result  has  been  obtained  by  Matthewson  (30)  for  thin 
compliant  coatings  with  Poisson's  ratios  in  the  range  0.2  to  0.3.  However, 
his  model  showed  that  the  inferfacial  shear  stress  increases  rapidly  for 
Poisson's  ratios  approaching  0.5. 

3.9  CONCLUSIONS 

It  has  been  shown  that  thin  hard  coatings  reduce  the  maximum  tensile 
stress  in  the  substrate  but  also  increase  the  maximum  stress  in  the  compo¬ 
site  structure  compared  to  the  values  for  an  uncoated  halfspace.  Thin 
hard  coatings  can  therefore  be  used  to  protect  brittle  materials  against 
elastic  contact  damage.  Although  the  stress  reduction  is  a  very  important 
indication  about  the  coating  performance,  more  specific  predictions  can 
only  be  made  when  information  is  available  on  other  factors  such  as  the 
fracture  toughness  of  the  coating  and  substrate  material,  the  interfacial 
strength,  pre-existing  stresses  and  the  flaw  size  distribution. 

With  regard  to  the  selection  of  optimum  coating  materials  the  present 
analysis  has  shown  that  materials  with  a  high  Young's  modulus  give  the 
largest  reduction  in  the  substrate  stresses  and  should  therefore  be  pre¬ 
ferred  provided  that  the  mechanical  properties  of  the  coating  material  are 
good  enough  to  withstand  the  high  stresses  which  develop  in  the  coating. 
Although  the  range  of  Poisson's  ratios  for  hard  materials  is  rather  small, 
materials  with  a  relatively  high  Poisson’s  ratio  (about  0.3)  should  be 
chosen  since  both  the  maximum  substrate  and  the  maximum  coating  stress 
decreases  with  increasing  Poisson's  ratio. 

When  considering  the  optimum  coating  thickness  it  should  be  noted 
that  the  analysis  indicates  two  regions  of  suitable  coating  thicknesses, 
namely  very  thin  (d/a  <  0,05)  or  relatively  thick  coatings  (d/a >0.4). 

For  the  relatively  thick  coatings  the  maximum  substrate  stresses  are  vir¬ 
tually  independent  of  the  coating  thickness  while  the  maximum  coating 
stresses  decrease  with  increasing  coating  thickness.  However,  for  these 
thick  coatings  the  problem  of  spontaneous  coating  debonding  due  to  residual 
stresses  is  a  serious  one.  For  the  thin  coatings  a  sizeable  reduction  in 
the  substrate  stresses  is  obtained  without  too  large  an  increase  in  the 
maximum  coating  stress.  Ball  indentation  experiments  on  carbon-coated 
germanium  by  van  der  Zwaag  and  Field  (43)  have  shown  that  it  is  indeed 
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possible  to  obtain  a  significant  increase  in  the  critical  load  for 
nucleation  using  thin  rigid  coatings. 


SECTION  IV 


THE  EFFECT  OF  THIN  HARD  COATINGS  FOR  A  LOCALLY  APPLIED  UNIFORM 
PRESSURE 

4.1  INTRODUCTION 

In  this  section  we  will  examine  the  effect  of  coating  thickness  and 
the  elastic  properties  of  the  coating  material  on  the  surface  and  inter¬ 
face  stresses  for  the  case  of  a  locally  applied  uniform  pressure.  The 
same  finite  element  model  as  used  in  the  previous  section  was  employed 
apart  from  a  change  in  the  prescribed  loading  condition  :  a  uniform  pre¬ 
ssure  over  the  chosen  contact  radius  instead  of  a  spherical  indentation. 

We  will  also  compare  the  static  and  transient  stress  field  in  a  un¬ 
iform  halfspace  due  to  such  a  loading  condition  to  find  the  limitations  of 
a  static  model  in  describing  impact  phenomena. 

4.2  RESULTS  FOR  AN  EXAMPLE  PROGRAM 

The  same  coating  parameters  as  used  in  the  example  program  in 
section  3.4  were  chosen  to  demonstrate  the  effect  of  a  hard  coating  on  the 
stress  field  due  to  a  local  uniform  pressure  :  E^  =  610  GPa;  =  0.20; 

Eg  =  61  GPa  and  Vg  =  0.25  and  the  ratio  of  the  coating  thickness  to  the 
contact  radius  =  0.05. 

The  calculated  radial  stress  distribution  at  the  coating  interface 
is  shown  in  Figure  22.  The  radial  stress  is  normalised  to  the  pressure 
prescribed  to  the  surface  nodes  within  the  contact  zone  area.  The  finite 
element  results  are  indicated  with  closed  circles.  The  dashed  line  in 
the  figure  indicates  the  radial  stress  in  a  uniform  halfspace  at  the  depth 
of  the  interface.  Even  for  this  small  coating  thickness  the  coating  has 
a  considerable  effect  on  the  radial  stress  field  both  within  and  outside 
the  contact  area.  The  same  applies  to  the  hoop  stress  field  which  is 
plotted  in  Figure  23  for  both  the  coated  and  uncoated  condition.  The  ef¬ 
fect  of  such  a  thin  coating  on  the  magnitude  of  the  axial  stresses  (not 
shown) is  however  negligible  for  the  present  loading  condition.  The  shear 
stress  field  at  the  interface  is  shown  in  Figure  24.  The  figure  shows  that 
there  is  a  large  shear  stress  which  reaches  its  maximum  at  the  contact  ra¬ 
dius.  The  radial  stress  at  the  coating  surface  is  plotted  in  Figure  25. 
For  this  small  coating  thickness  the  increase  in  the  radial  stress  at  the 
coating  surface  is  already  considerable  and  is  due  to  the  relatively  large 
deflections  of  the  substrate. 

In  conclusion,  the  picture  which  emerges  from  this  example  is  quali¬ 
tatively  the  same  as  for  the  example  program  in  the  previous  section. 
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4.3  EFFECT  OF  THE  COATING  PARAMETERS  ON  THE  RADIAL  STRESS 

The  radial  stress  distribution  at  the  interface  is  shown  in  Figure 
26  for  three  different  Young's  moduli  of  the  coating  material  (2,5,10  times 
that  of  the  substrate).  The  Poisson's  ratio  of  all  coatings  is  0.25  and 
the  coating  thickness  is  0.1  a.  The  same  tendency  is  observed  as  for  the 
ball  indentation  calculations  presented  in  the  previous  section;  namely, 
a  reduction  of  the  radial  stress  outside  the  contact  radius  with  increasing 
coating  stiffness  and  an  increase  in  the  radial  distance  at  which  the  max¬ 
imum  radial  stress  is  reached. 

The  effect  of  all  the  coating  parameters  on  the  maximum  radial  stress 
at  the  interface  between  coating  and  substrate  is  summarised  in  Figure  27. 

The  effect  of  the  elastic  properties  of  the  coating  and  its  thick¬ 
ness  on  the  maximum  tensile  stress  in  the  coating  surface  is  plotted  in 
Figure  28.  As  is  already  shown  in  Figure  25  this  maximum  occurs  at  or  near 
the  edge  of  the  loaded  area.  Also  for  the  present  loading  conditions  the 
maximum  tensile  stress  at  the  surface  is  considerably  higher  than  that  for 
the  uncoated  substrate.  When  comparing  the  thickness  dependence  of  the 
maximum  surface  stress  with  that  for  a  spherical  indentation  (Figure  19)  it 
will  be  clear  that  a  uniform  pressure  over  the  contact  area  causes  a  larger 
increase  in  the  surface  stresses  and  that  the  coating  thickness  at  which 
the  highest  tensile  stresses  occur  is  smaller.  These  differences  are  due 
to  the  differences  in  the  magnitude  of  the  pressure  near  the  edge  of  the 
contact  zone.  The  influence  of  the  elastic  properties  of  the  coating  mater¬ 
ial  on  these  curves  is  comparable  for  the  two  loading  conditions. 

4.4  DYNAMIC  AND  STATIC  STRESS  FIELD  DUE  TO  A  UNIFORM  PRESSURE 

The  solution  of  the  transient  stress  field  in  a  uniform  halfspace 
due  to  a  uniform  pressure  over  an  expanding  contact  area  has  been  determined 
by  Blowers  (44).  His  model  is  a  simplification  of  the  pressure  distribution 
in  a  solid  target  due  to  an  impacting  drop.  Recent  theoretical  (4,45,46) 
and  experimental  work  (47, 4")  has  shovTi  that  the  pressure  which  develops 
during  the  initial  stages  of  the  impact  is  not  very  uniform  but  rises 
sharply  towards  the  edge  of  the  contact  area.  However,  Adler  (49)  concludes 
from  the  good  agreement  between  the  Blowers  result  and  a  more  complex  finite 
difference  model  in  which  both  the  liquid  drop  and  the  target  are  considered 
(46),  that  the  contribution  of  the  edge  pressure  to  the  transient  radial 
stress  field  is  very  small. 

The  normalised  radial  stress  field  calculated  according  to  Blowers 
model  depends  on  both  the  material  properties  and  the  impact  velocity; 


unlike  the  static  case  where  the  normalised  stress  field  depends  on  the 
material  properties  only.  Examples  of  the  transient  radial  stress  field 
(normalised  to  the  impact  pressure  P)  are  shown  in  Figure  29  for  a  200  m  s 
and  a  600  m  s  *  drop  impact  on  zinc  sulphide  respectively.  It  will  be 
clear  from  these  figures  that  the  spatial  distribution  of  the  transient 
radial  stress  is  very  different  from  that  for  the  static  case  where  the 
surface  stress  outside  the  contact  area  is  given  by  (50) 
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Furthermore,  the  figures  show  that  the  ratio  of  the  maximum  tensile  stress 
and  the  applied  pressure  is  impact  velocity  dependent;  for  the  static  case 
the  ratio  is  0.25,  for  an  impact  velocity  of  200  ms',  0.9  and  for  a 
600  m  s  impact  1.6.  This  increase  is  due  to  the  fact  that  with  increas¬ 
ing  impact  velocity  the  contact  area  expands  faster.  Since  the  stress  wave 
velocity  in  the  material  remains  constant*  the  pulse  energy  becomes  more  lo¬ 
calised,  i.e.  the  resulting  stresses  are  increased. 

A. 5  DISCUSSION 

It  has  been  shown  that  for  the  present  loading  condition,  i.e.  a  lo¬ 
cally  applied  uniform  pressure,  hard  coatings  can  reduce  the  maximum  tensile 
stress  in  the  substrate  and  that  the  reduction  depends  primarily  on  the  coat¬ 
ing  thickness  and  the  Young's  modulus  of  the  coating  material. 

The  effect  of  the  coating  parameters  on  the  maximum  tensile  stress 
both  in  the  coating  and  the  substrate  is  about  the  same  as  for  the  spherical 
indentation;  in  particular  the  ranking  of  the  optimum  elastic  properties  of 
the  coating  material  is  identical. 

The  comparison  between  the  static  and  dynamic  stress  field  in  a  uni¬ 
form  halfspace  due  to  a  uniform  pressure  has  shown  that  there  are  consider¬ 
able  differences  both  in  the  distribution  and  the  magnitude  of  the  radial 
stresses.  The  differences  between  the  static  and  dynamic  stress  fields  will 
be  even  larger  for  layered  structures  since  in  the  second  case  the  final 
stress  field  will  be  the  sum  of  the  surface  or  Rayleigh  wave,  the  interface 
or  Stoneley  wave  and  the  various  incident  and  reflected  bulk  stress  waves. 
Furthermore,  the  type  of  stress  wave  will  depend  on  the  ratio  of  the  wave¬ 
length  and  the  coating  thickness.  Although  equations  for  various  elastic 
waves  in  layered  structures  have  been  developed  (see  for  a  general  intro¬ 
duction  Ewing  et  al  (51)) the  complete  solution  of  the  simultaneous  effect 
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of  all  the  stress  waves  on  the  final  stress  field  at  the  interface  and  in 
the  coating  will  require  considerable  effort;  in  particular  for  loading 
conditions  as  connlex  as  during  liquid  drop  impact. 

The  relevance  of  the  present  work  with  respect  to  protection 
against  rain  erosion  damage  should  therefore  be  seen  in  the  light  of  pro¬ 
tecting  the  substrate  from  contact  damage  during  handling  and  exposure  to 
(quasi-static)  dusty  atomospheres ,  since  a  reduction  of  the  initial  flaw 
size  will  cause  a  reduction  in  the  high  velocity  rain  impact  damage  in 
brittle  materials,  as  demonstrated  in  the  model  presented  in  Section  2. 
However,  recent  experiments  on  carbon-coated  germanium  (43)  have  shown 
that  hard  coatings  can  provide  protection  against  contact  damage  for 
static  (ball  indentations)  as  well  as  dynamic  (liquid  jet  impact)  loading 
conditions . 
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Failure  Probability 


Fracture  Stress  frequency 


lure  Probability 


Impact  Velocity 

5.  Theoretical  residual  strength  curves  for  soda-lime 
glass  (see  text), 

0  'as  received*  specimens' 

•  'abraded*  specimens. 


Impact  Velocity 

6.  Theoretical  residual  strength  curves  for  soda-lime 
glass  (  0  )  and  'tough-material'  (  •  ). 
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Finite  element  distribution  for  Hertzian  indentation 
on  a  coated  half space. 

Detail  of  7a.  showing  the  coating  elements  (shaded) and 
the  substrate  elements. 
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9.  Circumferential  stress  at  the  coating  interface 
(see  text) . 


Substrate  Halfspace 


Interfacial  shear  stress  at  the  interface  (see  text) 


Substrate  Halfspace 


Radial  Stress/  P*  .  Radial  Stress 


Figure 


16. 


Radial  stress  at  the  interface  for  four  different 
Young's  moduli  of  the  coating  material. 

(coating  thickness  ■>  0.1a  and  v  -  0.25). 
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Shear  Stress 


Figure 


20.  Variation  of  the  radial  shear  stress  as  a  function 
of  the  coating  thickness  (E  «>  10  E  and  v  ■>0.20 

O  8  C 

for  all  coatings). 

0  :  d  ■>  0.025a  •  :  d  ■■  0.05a 


Radial  stress  at  the  surface  for  a  uniform  pressure 
Dashed  line  indicates  the  radial  stress  for  an 
uncoated  substrate. 
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Figure  29.  Development  of  the  radial  stress  at  the  surface  of 
uncoated  zinc  sulphide  for  an  expanding  uniform 
pressure 

a)  impact  velocity  200  m  s  * 

b)  impact  velocity  600  m  s  * 
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